A systematic search for massive black hole binaries in SDSS spectroscopic sample

We present the results of a systematic search for massive black hole binaries in the Sloan Digital Sky Survey spectroscopic database. We focus on bound binaries, under the assumption that one of the black holes is active. In this framework, the broad lines associated to the accreting black hole are expected to show systematic velocity shifts with respect to the narrow lines, which trace the rest-frame of the galaxy. For a sample of 54586 quasars and 3929 galaxies at redshifts 0.1<z<1.5 we brute-force model each spectrum as a mixture of two quasars at two different redshifts. The spectral model is a data-driven dimensionality reduction of the SDSS quasar spectra based on a matrix factorization. We identified 32 objects with peculiar spectra. Nine of them can be interpreted as black hole binaries. This doubles the number of known black hole binary candidates. We also report on the discovery of a new class of extreme double-peaked emitters with exceptionally broad and faint Balmer lines. For all the interesting sources, we present detailed analysis of the spectra, and discuss possible interpretations.Comment: 10 pages, 2 figures, accepted for publication in Ap

Structure of the SMC - Stellar component distribution from 2MASS data

The spatial distribution of the SMC stellar component is investigated from 2MASS data. The morphology of the different age populations is presented. The center of the distribution is calculated and compared with previous estimations. The rotation of the stellar content and possible consequence of dark matter presence are discussed. The different stellar populations are identified through a CMD diagram of the 2MASS data. Isopleth contour maps are produced in every case, to reveal the spatial distribution. The derived density profiles are discussed. The older stellar population follows an exponential profile at projected diameters of about 5 kpc (~5 deg) for the major axis and ~4 kpc for the minor axis, centred at RA: 0h:51min, Dec: -73deg 7' (J2000.0). The centre coordinates are found the same for all the different age population maps and are in good accordance with the kinematical centre of the SMC. However they are found considerably different from the coordinates of the centre of the gas distribution. The fact that the older population found on an exponential disk, gives evidence that the stellar content is rotating, with a possible consequence of dark matter presence. The strong interactions between the MCs and the MilkyWay might explain the difference in the distributions of the stellar and gas components. The lack in the observed velocity element, that implies absence of rotation, and contradicts with the consequences of exponential profile of the stellar component, may also be a result of the gravitational interactions.Comment: 7 Pages, 6 figures, accepted for publication in A&

The nature of massive black hole binary candidates - II. Spectral energy distribution atlas

Recoiling supermassive black holes (SMBHs) are considered one plausible physical mechanism to explain high velocity shifts between narrow and broad emission lines sometimes observed in quasar spectra. If the sphere of influence of the recoiling SMBH is such that only the accretion disc is bound, the dusty torus would be left behind, hence the SED should then present distinctive features (i.e. a mid-infrared deficit). Here, we present results from fitting the spectral energy distributions (SEDs) of 32 type-1 AGN with high velocity shifts between broad and narrow lines. The aim is to find peculiar properties in the multiwavelength SEDs of such objects by comparing their physical parameters (torus and disc luminosity, intrinsic reddening, and size of the 12 μm emitter) with those estimated from a control sample of ∼1000 typical quasars selected from the Sloan Digital Sky Survey in the same redshift range. We find that all sources, with the possible exception of J1154+0134, analysed here present a significant amount of 12 μm emission. This is in contrast with a scenario of an SMBH displaced from the centre of the galaxy, as expected for an undergoing recoil event

Luminous AGB stars in nearby galaxies. A study using Virtual Observatory tools

Aims. This study focuses on very luminous Mbol<-6.0 mag AGB stars with J-Ks>1.5 mag and H-Ks>0.4 mag in the LMC, SMC, M31, and M33 from 2MASS data. Methods.The data were taken from the 2MASS All-Sky Point Source catalogue archive. We used Virtual Observatory tools and took advantage of its capabilities at various stages in the analysis. Results. It is well known that stars with the colors we selected correspond mainly to carbon stars. Although the most luminous AGBs detected here contain a large number of carbon stars,they are not included in existing catalogues produced from data in the optical domain, where they are not visible since they are dust-enshrouded. A comparison of the AGB stars detected with combined near and mid-infrared data from MSX and 2MASS in the LMC shows that 10% of the bright AGB stars are bright carbon stars never detected before and that the other 50% are OH/IR oxygen rich stars, whereas the 40% that remain were not cross-matched. Conclusions. The catalogues of the most luminous AGB stars compiled here are an important complement to existing data. In the LMC, these bright AGB stars are centrally located, whereas they are concentrated in an active star-formation ring in M31. In the SMC and M33, there are not enough of them to draw definite conclusions, although they tend to be centrally located. Their luminosity functions are similar for the four galaxies we studied.Comment: 16 pages, 12 figures, 4 tables (Appendix A), accepted in A&

The expected performance of stellar parametrization with Gaia spectrophotometry

Gaia will obtain astrometry and spectrophotometry for essentially all sources in the sky down to a broad band magnitude limit of G=20, an expected yield of 10^9 stars. Its main scientific objective is to reveal the formation and evolution of our Galaxy through chemo-dynamical analysis. In addition to inferring positions, parallaxes and proper motions from the astrometry, we must also infer the astrophysical parameters of the stars from the spectrophotometry, the BP/RP spectrum. Here we investigate the performance of three different algorithms (SVM, ILIUM, Aeneas) for estimating the effective temperature, line-of-sight interstellar extinction, metallicity and surface gravity of A-M stars over a wide range of these parameters and over the full magnitude range Gaia will observe (G=6-20mag). One of the algorithms, Aeneas, infers the posterior probability density function over all parameters, and can optionally take into account the parallax and the Hertzsprung-Russell diagram to improve the estimates. For all algorithms the accuracy of estimation depends on G and on the value of the parameters themselves, so a broad summary of performance is only approximate. For stars at G=15 with less than two magnitudes extinction, we expect to be able to estimate Teff to within 1%, logg to 0.1-0.2dex, and [Fe/H] (for FGKM stars) to 0.1-0.2dex, just using the BP/RP spectrum (mean absolute error statistics are quoted). Performance degrades at larger extinctions, but not always by a large amount. Extinction can be estimated to an accuracy of 0.05-0.2mag for stars across the full parameter range with a priori unknown extinction between 0 and 10mag. Performance degrades at fainter magnitudes, but even at G=19 we can estimate logg to better than 0.2dex for all spectral types, and [Fe/H] to within 0.35dex for FGKM stars, for extinctions below 1mag.Comment: MNRAS, in press. Minor corrections made in v

Pulsar timing arrays and the challenge of massive black hole binary astrophysics

Pulsar timing arrays (PTAs) are designed to detect gravitational waves (GWs) at nHz frequencies. The expected dominant signal is given by the superposition of all waves emitted by the cosmological population of supermassive black hole (SMBH) binaries. Such superposition creates an incoherent stochastic background, on top of which particularly bright or nearby sources might be individually resolved. In this contribution I describe the properties of the expected GW signal, highlighting its dependence on the overall binary population, the relation between SMBHs and their hosts, and their coupling with the stellar and gaseous environment. I describe the status of current PTA efforts, and prospect of future detection and SMBH binary astrophysics.Comment: 18 pages, 4 figures. To appear in the Proceedings of the 2014 Sant Cugat Forum on Astrophysics. Astrophysics and Space Science Proceedings, ed. C.Sopuerta (Berlin: Springer-Verlag

Fitting the integrated Spectral Energy Distributions of Galaxies

Fitting the spectral energy distributions (SEDs) of galaxies is an almost universally used technique that has matured significantly in the last decade. Model predictions and fitting procedures have improved significantly over this time, attempting to keep up with the vastly increased volume and quality of available data. We review here the field of SED fitting, describing the modelling of ultraviolet to infrared galaxy SEDs, the creation of multiwavelength data sets, and the methods used to fit model SEDs to observed galaxy data sets. We touch upon the achievements and challenges in the major ingredients of SED fitting, with a special emphasis on describing the interplay between the quality of the available data, the quality of the available models, and the best fitting technique to use in order to obtain a realistic measurement as well as realistic uncertainties. We conclude that SED fitting can be used effectively to derive a range of physical properties of galaxies, such as redshift, stellar masses, star formation rates, dust masses, and metallicities, with care taken not to over-interpret the available data. Yet there still exist many issues such as estimating the age of the oldest stars in a galaxy, finer details ofdust properties and dust-star geometry, and the influences of poorly understood, luminous stellar types and phases. The challenge for the coming years will be to improve both the models and the observational data sets to resolve these uncertainties. The present review will be made available on an interactive, moderated web page (sedfitting.org), where the community can access and change the text. The intention is to expand the text and keep it up to date over the coming years.Comment: 54 pages, 26 figures, Accepted for publication in Astrophysics & Space Scienc

Supermassive Black Hole Binaries: The Search Continues

Gravitationally bound supermassive black hole binaries (SBHBs) are thought to be a natural product of galactic mergers and growth of the large scale structure in the universe. They however remain observationally elusive, thus raising a question about characteristic observational signatures associated with these systems. In this conference proceeding I discuss current theoretical understanding and latest advances and prospects in observational searches for SBHBs.Comment: 17 pages, 4 figures. To appear in the Proceedings of 2014 Sant Cugat Forum on Astrophysics. Astrophysics and Space Science Proceedings, ed. C.Sopuerta (Berlin: Springer-Verlag

<i>Gaia</i> Data Release 1. Summary of the astrometric, photometric, and survey properties

Context. At about 1000 days after the launch of Gaia we present the first Gaia data release, Gaia DR1, consisting of astrometry and photometry for over 1 billion sources brighter than magnitude 20.7. Aims. A summary of Gaia DR1 is presented along with illustrations of the scientific quality of the data, followed by a discussion of the limitations due to the preliminary nature of this release. Methods. The raw data collected by Gaia during the first 14 months of the mission have been processed by the Gaia Data Processing and Analysis Consortium (DPAC) and turned into an astrometric and photometric catalogue. Results. Gaia DR1 consists of three components: a primary astrometric data set which contains the positions, parallaxes, and mean proper motions for about 2 million of the brightest stars in common with the HIPPARCOS and Tycho-2 catalogues – a realisation of the Tycho-Gaia Astrometric Solution (TGAS) – and a secondary astrometric data set containing the positions for an additional 1.1 billion sources. The second component is the photometric data set, consisting of mean G-band magnitudes for all sources. The G-band light curves and the characteristics of ∼3000 Cepheid and RR-Lyrae stars, observed at high cadence around the south ecliptic pole, form the third component. For the primary astrometric data set the typical uncertainty is about 0.3 mas for the positions and parallaxes, and about 1 mas yr−1 for the proper motions. A systematic component of ∼0.3 mas should be added to the parallax uncertainties. For the subset of ∼94 000 HIPPARCOS stars in the primary data set, the proper motions are much more precise at about 0.06 mas yr−1. For the secondary astrometric data set, the typical uncertainty of the positions is ∼10 mas. The median uncertainties on the mean G-band magnitudes range from the mmag level to ∼0.03 mag over the magnitude range 5 to 20.7. Conclusions. Gaia DR1 is an important milestone ahead of the next Gaia data release, which will feature five-parameter astrometry for all sources. Extensive validation shows that Gaia DR1 represents a major advance in the mapping of the heavens and the availability of basic stellar data that underpin observational astrophysics. Nevertheless, the very preliminary nature of this first Gaia data release does lead to a number of important limitations to the data quality which should be carefully considered before drawing conclusions from the data